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, Abstract. Inclusive spectra of charged particles at midrapidity in Au+Au 

collisions at y/s^ = 130 GeV and 200 GeV were measured with the STAR 
t-H , detector at RHIC. The measured mean transverse momentum {pr) shows a 

characteristic dependence on charged particle multiplicity and beam energy in 
Au+Au collisions that is distinctly different from pp, pp and e + e~ collisions. 
A 32% ± 3%(syst) increase in (pr) from pp to Au+Au collisions was observed 
ly-j , at 200 GeV. While the charged multiplicity was found to increase by 19% ± 

• 5%(syst) from ^/s^ — 130 GeV to 200 GeV, no significant difference in (pr) 

C*~) , was found between the two energies. A comparison with model predictions is 

discussed. 
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j> \ 1. Introduction 

•i-H . 

' The systematic study of charged particle multiplicity in Au+Au collisions at center 

of mass energies y/s^— 130 GeV and 200 GeV was conducted by all RHIC ex- 
periments in great detail [Ql El EJ El EJ- These studies were accompanied by 
intensive theoretical work successfully connecting the produced particle densities 
with the initial conditions in these collisions. Essentially all models that attempt to 
describe the evolution of matter produced in relativistic heavy-ion collisions predict 
a strong correlation between the average transverse momentum and the multiplicity 
of the event, typically (pr) oc dN^/dr). This simply reflects the basic thermody- 
namic concept inherent to the models: large i\T c h means large entropy, large entropy 
means higher temperature T, which, following the general concepts of Landau's 
hydrodynamical model, implies an increase in (pr)- The use of thermodynamics 
in this argument does not, of course, preclude a dynamical description in terms of 
(non-)perturbative Quantum Chromodynamics (QCD). It simply takes advantage 
of the fact that if the actual physical states of the system are close to a canoni- 
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cal distribution, thermodynamics is an ideal tool to analyze the measurements and 
relate them to bulk properties of hot matter. 

One prominent class of models that allows a direct coupling of global observ- 
ables with the initial conditions of the collision are those based on the concept of 
gluon saturation. In these models particle production is related to the initial gluon 
distribution. A semi-classic calculation can be carried out due to the high gluon 
density, a saturation phenomenon also known as the Color Glass Condensate [0 EJ- 
The only scale that controls the gluon distribution at a given Bjorken x is the sat- 
uration momentum scale Q s . If the gluons realize themselves directly to hadrons, 
the transverse momentum of the hadrons will be directly related to the saturation 
scale Q s as well as the multiplicity of the particle production[|21 [5]. In fact, this 
assumption was used in some predictions that successfully predicted the particle 
multiplicity along with their centrality and pseudorapidity dependence [[3J. 

So far, however, only few studies were undertaken to verify if these models are 
also able to describe the characteristics of (pr) at RHIC. In this article we present 
results on the characteristics of (i) the charged particle densities and (ii) of (pr) 
using data recorded by the STAR experiment. 

2. Inclusive charged hadron spectra 

The data reported here were taken with the Solenoidal Tracker at RHIC (STAR) [ 
I10j . The main tracking device in STAR is a large Time Projection Chamber (TPC) 
which provides momentum information and particle identification for charged parti- 
cles by measuring their ionization energy loss. A Central Trigger Barrel (CTB) con- 
structed of scintillator paddles surrounding the TPC and two Zero Degree Calorime- 
ters (ZDC) were used for triggering. Data were taken for Au+Au collisions at 
^^"=130 GeV in the summer of 2000 (run I) and y /s^=200 GeV in the fall of 
2001 (run II). The field setting for run I was at 0.25 T while at run II it was at 
the nominal field setting of 0.5 T. The increased field resulted in a factor of three 
improvement in momentum resolution. In both runs a minimum bias trigger was 
defined by using coincidences between the two ZDCs located ±18 meters away from 
the interaction point. For details see Ref. [Illj. 

To measure the differential invariant cross-section of the produced charged par- 
ticle for different centralities tracks with transverse momenta between 0.2 and 2.0 
GeV/c were selected. The reconstruction efficiency was determined by embedding 
simulated tracks into real events at the raw data level, reconstructing the full events, 
and comparing the simulated input to the reconstructed output. Charged hadrons 
(h + + h _ )/2 were approximated by the summed yields of , K^, p and p. The 
efficiencies for the identified particles are between 70% and 80% for central events 
at pr > 0.5 GeV/c. At lower pr (0.2 to 0.5 GeV/c), the efficiency for pion is 
about 70% and the efficiencies of kaon and antiproton are from 30% to about 70% 
due to energy loss, absorption and decay in the detector. The final efficiency of the 
charged particle spectrum is the weighted efficiency of all tt , K ± and p ± according 
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Fig. 1. Inclusive invariant pr distributions of (h + + h~)/2 for ^/s ^ 
GeV Au+Au collisions. The lines are power law fits to the data (see text). 

to their yield ratios measured at ^/s NN — 130 GeV and the proper Jacobian factor 
dy/dij(pT 7 rj). It is almost identical to that of the pions for all px- 

Background particles mis-identified as primary particles originate from various 
sources. The dominant ones are: hadrons from weak decays of Kg, A and A, weak 
decays of primary charged particles, and hadrons from secondary interactions in 
the detectors. The background corrections for the former were calculated by using 
measured yields and spectra. The later two are important only at low px and 
we used HIJING with a detailed detector simulation to reproduce the background 
spectra for each centrality. 

Figure n shows the fully corrected inclusive pr distributions of primary charged 
particles for y / s^=130 and 200 GeV Au+Au collisions for 10 and 7 centrality 
bins, respectively. The spectra were fit with a power law function of the form 
A(l + pr/po) 71 - The average momentum is defined as (px) — 2 po/(n ~ 3). To 
extract the (jot) we substitute po by (pr)(n — 3)/2 and fit the spectrum with (jot) 
and n as a free parameters. The difference between the integral of the fit function 
and the sum of the data points in the fiducial kinematic range is less than 1%. 
The extrapolation to the total pt is 21-26% from central to peripheral centrality. 
The systematical uncertainty of (pt) was estimated to 2%. The function, fit in 
the range of 0.2 GeV/c to 2.0 GeV/c, was then extrapolated to the full transverse 
momentum range to get the total iVch- The systematical uncertainty on N c h is 7%. 

For the 5% most central collisions. The charged hadron multiplicity at ^/s NN =200 
GeV (130 GeV) is di\W d »7l|>?|<0.5 = 690 ± 55(syst) (579 ± 60(syst)) and the 
mean transverse momentum (p T ) is 0.517 ± 0.012GeV/c (0.519 ± O.OlGeV/c). For 
^/s NN =200 GeV the multiplicity per participant pair increases by 38% relative to 
pp and 52% compared to nuclear collisions at v / s^=17 GeV. The (pt) increases 
by 32% and 20% respectively from 0.392 GeV/c in pp collisions and ~ 0.429 GeV/c 
in Pb+Pb collisions at SPS [ITT1 IT2]. 
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3. Centrality dependence of charged particle densities 

The RHIC experiments all have different methods of determining centrality. In 
general, all measure some final state variable (or combination of these variables) 
that can be related to fractions of the total measured hadronic cross-section. These 
centrality bins are then mapped to variables that allow for direct comparison be- 
tween experiments. The primary variables that are used for comparison are (i) the 
number of participating nucleons in the collision A/p ar t and (ii) the number of bi- 
nary collisions Nhin, which are determined via the application of the Glauber model 
phenomenology. The basic concept of the Glauber model is to treat a Au+Au col- 
lision as a superposition of many independent nucleon-nucleon (N+N) collisions. 
Thus, the only parameters that the model depends on are a nuclear density profile 
(Woods-Saxon) and the non-diffractive inelastic N+N cross-section. The former is 
well measured in e+Au scattering experiments. The latter is well established by 
many previous experiments. With this, it is possible to determine both the total 
number of independent N+N collisions (A^in) and the total number of nucleons 
that participate in the collisions (N part ) as a function of impact parameter (b). 
There are two separate implementations of the Glauber model: Optical and Monte 
Carlo. In the optical formalism Abin and iVp ar t are directly determined by an an- 
alytic integration of overlapping Woods-Saxon distributions. In the M.C. method 
an arbitrary number of events are simulated using a computer program and the 
resulting distributions of AT bin and N palt are determined. Once the distributions 
dc/d/Vpart and da/dNun are determined, these histograms are binned according to 
fractions of the total cross-section. This determines the mean values of A^in and 
A^part for each centrality class. In this context it should be noted that in the M.C. 
approach one requires that all nucleons in either nucleus be separated by a distance 
d > rf m in where d m - ln = 0.4 fm is characteristic of the length of the repulsive nucleon- 
nucleon force. Consequently the total Au+Au cross section predicted by the Monte 
Carlo approach depends strongly on d m i n - A value of c? m ; n = 0.4 fm yields a total 
cross-section of 7.2 barns, while d min = fm yields 6.8 barns. 

While at SPS energies both approaches yielded consistent answers, at RHIC 
energies, the optical approach may yield systematically smaller values of the total 
Au+Au cross-section. Therefore, any results reported in terms of Glauber quantities 
must be carefully interpreted based upon specifics of the Glauber calculations. 

In the case of STAR the centrality was determined by the multiplicity of all 
charged particles within \r/\ < 0.5 and then expressed in terms of the fractional total 
hadronic cross section. The measured data, i.e. the dN/dN c h distributions, are 
then mapped to the corresponding distribution obtained from Glauber calculations 
as described above, thus relating Ap art and Abin to the measured distributions. 

Figure El shows the extracted yields as the pseudo-rapidity density per par- 
ticipant pair 2/Afp ar t dN/d-q versus the number of participants Ap ar t for Au+Au 
collisions at y's NN = 130 GeV, where for the left plot we used a Monte Carlo Glauber 
calculation to derive AT part while for the right plot we used an optical Glauber model. 

The parameterization of pseudo-rapidity density of the EKRT model [EI1> a 
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Fig. 2. Rapidity density per participant nucleon pair 2/(iV part ) dN/dij versus 
number of participants (N par t). On the left plot iVp ar t was calculated using a MC 
Glauber approach while on the right plot an optical Glauber model was used. 



model based on the assumption of final state gluon saturation, fails to describe the 
data when using the Monte Carlo Glauber calculations, but agrees nicely with when 
using the optical Glauber calculation. Both, our data and the EKRT parameter- 
ization do not approach the charged particle pseudo-rapidity density n pp = 2.25 
(measured from non-single diffractive pp interactions and parameterized as n pp = 
2.5 ± 1.0 - (0.25 ± 0.19) ln(s) + (0.023 ± 0.008) ln 2 (s) [Ej) in the limit of N paIt =2. 
One might argue that the model breaks down for very peripheral collisions. The 
centrality dependence from mid-central to central collisions is in fact too weak to 
rule out the EKRT model. 

Comparing the two-component model described in Ref. [ I15j with our data 
shows the exact opposite behavior. Fitting the proposed function to our data we 
find agreement between the fit and our measurement only when using the Monte 
Carlo Glauber calculation (left plot in Fig.|2J). Here, we evaluate the fraction of hard 
collisions Xhard as 0.08 ±0.01 for our most central events. The fraction of produced 
particles originating from hard collisions is then calculated as Fhard = x n pp -/Vbm/ 
(dNch/di]) = 0.36 ±0.2. Both values Zhard and Fh ar d agree nicely with the numbers 
quoted by Ref. [0| and [ES|. 

It should be noted that in general the M.C. Glauber method is considered to 
be more accurate than the calculations in the optical limit. Nevertheless, in both 
approaches the systematic uncertainties in the calculation of iV part are considerable 
large for very peripheral collisions, the region where the differences between the 
various models become more apparent. For the centrality fraction of 70-80% the 
uncertainty on iV part already reaches 30%. 
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4. Characteristics of {pj) in heavy-ion collisions 

Using the spectra from Fig. ^ and the power law fits described above we now turn 
to the centrality dependence of (pr) ■ 

Figure 3 shows the (pt) as a func- 
tion of charge multiplicity for both en- 
ergies. The open symbols are from 200 
GeV data and the solid stars are mea- 
surements from 130 GeV data. The 
solid triangle depicts the UA1 measure- 
ment [ Hi . We also show the calcu- 
lation from Eq. ^ discussed in more 
details below. The simulations from 
RQMD with and without rescattering 
are shown as dotted and dashed lines, 
respectively. The dot-dashed line de- 
picts the simulation from HIJING with 
default settings. The (pr) increases as 
function of centrality from peripheral 
to semi-central collisions and seems to Fig. 3. (pr) as a function of N c h at y/s^= 
reach a plateau for central collisions. 130 GeV and 200 GeV. The solid line is a 

Figure 0] has two panels. The solid prediction assuming scaling of (px) with 7V c h 
triangles in the top panel depict the in- as described in the text. Also shown are the 
elusive charged particle ratios between simulations from RQMD with and without 
130 GeV and 200 GeV as a function of rescattering and HIJING 1.35. All model 
number of participants (neglecting the curves are for y / s^"= 130 GeV only, 
small difference of V part between these 

two beam energies). The open diamonds are measurements from the PHOBOS 
Collaboration For central collisions the charged particle multiplicity increases 
from 130 GeV to 200 GeV by 19%. This result is comparable with the results from 
other RHIC experiments However, the systematical uncertainty at peripheral 
collisions is too large to be conclusive as shown in the top panel. Surprisingly, one 
finds no increase of (pr) within the current systematic errors for any centrality bin 
as depicted in the bottom panel. 

Figure El shows (pt) of negatively charged particles h~ from NA49 and charged 
particles (h + + h~)/2 from various other experiments as function of ^/s NN ior pp, pp 
and central AA collisions. The dotted line is a parametrization of the measured (px) 
in pp collisions For the e + e~ data, the dot-dashed line is a prediction from 

JETSET [El- The fact that the (p T ) from AA collisions is distinctly different from 
both pp and e + e~ indicates that AA collisions are not simple superpositions of the 
elementary collisions. 

As already discussed in the introduction, saturation [[§] and hydrodynamical [ 
IT%] models predict some scaling behavior of (pr) 2 ~ (^)aa/ 7T ^ 2 - 

In order to get a quantitative analysis on the dependence between (pt) and 
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multiplicity, we require that the scaling dependence is linear and it should also 
satisfy the pp results. Thereby, there is only one free parameter as shown in Eq.^ 



{Pi 



" - * / ((pt)pp - a); 

Spp 



(1) 



where s AA — {dN / drf) AA / it R 2 and s pp = (dN /dij) pp /irr 2 are the multiplicity densi- 
ties per unit pseudorapidity per unit transverse area in AA and pp collisions (i? = 
tqA 1 ^ 3 ) and a — 0.3 GeV/c is a constant chosen to describe the AA data. For 
non-central collisions, the calculation of the transverse area is not straight forward. 
We take the centrality dependence parametrization of s AA from Ref. [EJ. The solid 
curves in Fig. 3 and Fig. [3] show the energy dependence and centrality dependence 
of (px) from Eq. The curves describe the data well from SPS to RHIC at 130 
GeV. However, the data show that the increase in (px) from ^/s NN = 130 to 200 GeV 
is not as strong as the thick-line indicates which may invalidate the scaling law from 
the saturation model . On the other hand, the energy dependence might indicate 
the importance of early thermalization in the partonic stage, as proposed in ref [ 
119) . Such initial partonic activity is consistent with the early development of flow 
as indicated from (V2) measurements at RHIC [ [2D] - It will be interesting to see if 
the description holds with an energy scan of lower beam energies at RHIC. 



5. Summary 

In summary, STAR has measured charged particle multiplicities and transverse 
momentum spectra at mid rapidity in Au+Au collisions at ^/s NN =130 GeV and 
200 GeV. Their centrality and beam energy dependences are distinctly different 
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from elementary collisions (e + e~, pp and pp). The saturation model with direct 
hadron production from gluons seems unable to explain the lack of increase in (jpr) 
from 130 GeV to 200 GeV. However, the saturation model with thermalization [ 
I19| might be consistent with our result. Further measurements with various beam 
energies and beam species such as d+Au may be able to further distinguish different 
production mechanisms. 
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